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ABSTRACT 

An  error  probability  analysis  is  performed  for  an  M-ary  orthogonal  frequency 
shift  keying  (MFSK)  receiver  employing  fast  frequency-hopped  (FH)  spread  spectrum 
waveforms  transmitted  over  a  fading  channel  with  partial-band  interference.  The 
partial-band  interference  is  modeled  as  a  Gaussian  process.  Wideband  thermal  noise  is 
also  included  in  the  analysis.  Diversity  is  performed  using  multiple  hops  per  data  bit. 
Each  diversity  reception  is  assumed  to  fade  independently  according  to  a  nonselective 
Rician  process.  A  nonlinear  combination  procedure  referred  to  as  self-normalization 
combining  is  employed  by  the  receiver  to  minimize  partial-band  interference  effects. 
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I.  INTRODUCTION 

Spread  spectrum  communication  systems  are  designed  to  use  signals  having  a 
much  wider  bandwidth  than  ordinary  communication  systems.  Some  of  advantages 
of  spread  spectrum  signals  are  the  energy  gain  achivable  against  a  narrow  band  jam- 
mer, low  probability  of  detection  (LPD),  low  probability  of  intercept  (LPI),  multiple 
access  operation,  and  the  capability  to  overcome  the  effects  of  multipath  [Ref.  1]. 
The  primary  disadvantages  of  spread  spectrum  signals  are  the  time  and  frequency 
synchronization  required  by  the  receiver  which  make  it  difficult  to  implement  the 
system  and  the  wide  frequency  band  requirement.  There  are  three  primary  spread 
spectrum  communication  methods:  time  hopping  (TH),  frequency  hopping  (FH)  and 
direct  spreading  (DS).  This  thesis  presents  an  error  probability  analysis  of  two  fast 
frequency-hopped  M-ary  orthogonal  frequency-shift  keyed  (FFH/MFSK)  systems, 
where  the  term  'fast'  implies  that  the  frequency-hopping  rate  is  greater  than  the 
symbol  rate.  The  first  to  be  examined  is  a  conventional  FFH/MFSK  system  and, 
the  second  is  a  self-normalization  FFH/MFSK  system.  Both  systems  use  noncoher- 
ent detection.  The  channel  is  assumed  to  experience  both  fading  and  partial-band 
interference.  The  performance  of  the  two  receivers  is  compared  with  each  other  as 
well  as  with  another  type  of  FFH/MFSK  receiver  referred  to  as  a  noise-normalization 
receiver  [Ref.  2]. 

A.      Conventional  Receiver  and  Self-Normalization  Receiver 

In  the  conventional  FFH/MFSK  receiver,  illustrated  in  Fig.  1.1,  the  input 
signal  passes  through  a  wideband  radio  frequency  (RF)  filter  and  is  mixed  with  the 
frequencies  which  are  created  pseudo-randomly  in  the  frequency  synthesizer  to  dehop 
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Figure  1.1:  Block  Diagram  of  Conventional  FFH/MFSK  noncoherent  Re- 
ceiver 


and  convert  the  signal  frequency  down  to  intermediate  frequency  (IF).  The  bandwidth 
of  the  RF  filter  (W)  must  be  wide  enough  to  pass  all  possible  hopping  frequencies, 
while  the  bandwidth  of  the  IF  filter  (B)  needs  pass  only  one  hopping  frequency. 
Thus,  the  bandwidth  of  the  RF  filter  (W)  is  approximately  equal  to  NB,  where  TV  is 
number  of  the  hopping  frequencies.  The  receiver  has  M  different  branches  for  each  of 
the  M  different  symbols.  The  IF  signal  is  demodulated  with  a  bandpass  filter  and  a 
quadratic  detector.  The  output  signal  from  each  quadratic  detector  is  sampled  L  times 
and  combined  to  provide  a  decision  statistic  between  the  M  available  outputs.  The 
transmitter  uses  fast  frequency-hopping  to  provide  L  independent  samples  through 
frequency  diversity.  The  receiver  is  assumed  to  dehop  the  signal  with  perfect  timing. 
In  the  self-normalization  receiver,  illustrated  in  Fig.  1.2,  each  of  the  M  quadratic 
detector  outputs  are  divided  by  the  sum  of  the  quadratic  detector  outputs  from  all 
M  channels.  This  normalization  factor  is  the  only  difference  from  the  conventional 
receiver.  When  partial-band  interference  present,  the  hops  that  are  corrupted  by  the 
partial-band  noise  jamming  dominate  the  performance  of  the  conventional  detector 
since  each  hop  is  equally  weighted  in  the  overall  decision  statistic.  The  normalization 
factor  minimizes  this  effect  in  the  self-normalization  receiver.  As  with  the  conven- 
tional FFH/MFSK  receiver,  the  L  samples  at  the  output  of  the  M  quadratic  detectors 
are  combined  to  provide  a  decision  statistic  between  the  M  available  quadratic  detec- 
tor outputs,  and  the  receiver  is  assumed  to  dehop  the  signal  with  perfect  timing. 

B.      Fading 

For  both  receivers  investigated  in  this  thesis,  the  received  signal  is  assumed 
to  be  a  combination  of  a  single,  nonfaded  (direct)  component  and  many  reflected 
(diffuse)  signals.  By  modelling  the  fading  as  Rician,  in  which  the  signal  is  considered 
to  consist  of  a  direct  signal  component  and  a  diffuse  signal  component,  we  obtain  a 


Figure  1.2:  Block  Diagram  of  Self-Normalization  FFH/MFSK  noncoherent 
Receiver 


general  result  that  is  valid  in  the  limit  of  large  direct-to-diffuse  signal  power  ratios 
for  Gaussian  channels  and  in  the  limit  of  small  direct-to-diffuse  signal  power  ratios 
for  Rayleigh  fading  channels  as  well  as  the  general  case  where  the  effects  of  both  the 
direct  and  diffuse  components  of  the  signal  must  be  included  in  the  analysis.  For  a 
Rician  faded  signal  with  amplitude  Vza^-,  the  probability  density  function  is  [Ref.  3] 

fAk(ak)  =  —  exp  [—^T-j  'o  {—)  u(ak)  (1.1) 

where  a\  is  the  power  of  the  direct  signal  component,  2<j2  is  the  power  of  the  diffuse 
signal  component,  «(•)  is  the  unit  step  function,  and  Iq{*)  is  the  zeroth-order  modified 
Bessel  function  of  the  first  kind.  When  the  direct  signal  component  power  is  zero, 
which  occurs  in  times  of  deep  channel  fading  on  the  channel  such  as  when  the  direct 
communication  path  is  blocked  by  terrain  or  other  obstacles,  the  probability  density 
function  reverts  to  the  Rayleigh  probability  density  function  [Ref.  3] 

fAM  =  ^exp(-^ju(ak)  (1.2) 

Also  in  this  thesis,  we  assume  that  the  smallest  spacing  between  frequency  hop 
slots  is  larger  than  the  coherence  bandwidth  of  the  channel  [Ref.  4].  As  a  result,  each 
hop  of  a  symbol  experiences  independent  fading.  In  addition,  the  channel  for  each 
hop  of  a  symbol  is  modeled  as  a  frequency-nonselective,  slowly  fading  Rician  process. 
This  implies  that  the  signal  bandwidth  is  much  smaller  than  the  coherence  bandwidth 
of  the  channel  and  that  the  hop  duration  is  much  smaller  than  the  coherence  time  of 
the  channel  [Ref.  5].  The  latter  assumption  is  equivalent  to  requiring  the  hop  rate 
to  be  large  compared  to  the  Doppler  spread  of  the  channel.  As  a  result,  the  signal 
amplitude  can  be  modeled  as  a  Rician  random  variable  that  remains  fixed  at  least 
for  the  duration  of  a  single  hop. 


C.      Partial-Band  Interference 

One  of  the  applications  of  frequency-hopped  spread  spectrum  signals  is  to  reduce 
receiver  performance  degradation  due  to  narrowband  interference  (either  intentional 
or  otherwise).  This  is  accomplished  by  transmitting  each  symbol  at  a  different  carrier 
frequency  according  to  some  apparently  random  pattern  known  only  by  the  trans- 
mitter and  the  receiver.  The  bandwidth  of  the  frequency-hopped  signal  is  generally 
much  larger  than  the  signal  bandwidth  in  the  absence  of  frequency-hopping.  Conse- 
quently, in  order  to  interfere  with  the  entire  spread  spectrum  band,  given  a  fixed  total 
interference  power  the  noise  power  spectral  density  of  the  narrowband  interference  is 
reduced.  Rather  than  interfere  with  the  entire  spread-spectrum  bandwidth  simultane- 
ously with  the  consequent  reduction  in  noise  power  spectral  density,  the  narrowband 
interference  may  affect  only  a  portion  of  the  total  spread  spectrum  bandwidth  at 
any  one  time,  and  the  portion  of  the  total  spread  spectrum  bandwidth  experiencing 
narrowband  interference  may  also  change  in  an  apparently  random  manner.  This 
is  referred  to  as  partial-band  interference.  In  the  case  of  partial-band  interference, 
only  some  of  the  transmitted  symbols  will  be  affected  by  narrowband  interference, 
and  the  question  arises  as  to  whether  performance  can  be  improved  by  implementing 
diversity  in  the  form  of  fast  frequency-hopping;  that  is,  each  information  symbol  is 
transmitted  multiple  times,  and  each  transmission  is  at  a  different  apparently  random 
carrier  frequency  within  the  overall  spread  spectrum  bandwidth. 

As  mentioned  previously,  in  this  thesis  the  performance  of  two  different  fast 
frequency-hopped  M-ary  orthogonal  frequency-shift  keyed  (FFH/MFSK)  noncoher- 
ent receivers  is  investigated.  The  performance  of  the  conventional  FFH/MFSK  non- 
coherent receiver  in  a  Rician  fading  channel  with  only  wideband  noise  (no  partial- 
band  interference)  has  been  previously  investigated  [Ref.  6].  For  the  self-normalized 
FFH/MFSK  receiver,  system  performance  is  obtained  as  an  upper  bound  on  the 
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probability  of  bit  error,  while  for  the  conventional  FFH/MFSK  receiver,  the  actual 
probability  of  bit  error  is  obtained.  For  both  systems,  the  transmitter  is  assumed  to 
implement 

In  addition  to  wideband  Gaussian  noise,  the  channel  is  assumed  to  be  affected 
by  narrowband  noise  in  the  form  of  partial-band  interference.  The  interference  that 
is  considered  in  this  thesis  may  be  due  to  a  partial-band  jammer  as  well  as  other 
unintended  narrowband  interferences.  The  partial-band  interference  is  modeled  as 
additive  Gaussian  noise  and  is  assumed  to  corrupt  only  a  fraction  7,  where  1  >  7  >  0, 
of  the  entire  spread  spectrum  bandwidth  at  any  one  time.  This  is  illustrated  in  Fig. 
1.3  where  5  represents  the  signal  power  level,  Nj/2  the  jamming  noise  plus  thermal 
noise  power  spectral  density,  iVo/2  the  thermal  noise  power  spectral  density,  B  the 
bandwidth  of  one  frequency  cell,  and  W  the  entire  frequency  band.  In  addition, 
partial-band  interference,  when  present,  is  assumed  to  be  present  in  each  branch  of  the 
MFSK  demodulator,  and  the  fraction  of  the  spread  spectrum  bandwidth  experiencing 
partial-band  interference  is  assumed  to  be  the  same  for  all  hops  of  a  symbol.  Hence, 
7  is  the  probability  that  partial-band  interference  is  present  in  all  M  channels  of  the 
receiver,  and  1  —  7  is  the  probability  that  partial-band  interference  is  not  present  in 
all  M  channels  of  the  receiver.  If  Nj/2  is  defined  to  be  the  average  power  spectral 
density  of  narrowband  interference  over  the  entire  spread  spectrum  bandwidth,  then 
7_1  Nj/2  is  the  power  spectral  density  of  partial-band  interference  when  it  is  present. 
Thermal  noise  and  other  wideband  noise  that  corrupt  the  channel  are  modeled  as 
additive  white  Gaussian  noise,  and  the  power  spectral  density  of  this  wideband  noise  is 
defined  as  N0/2.  Thus,  the  power  spectral  density  of  the  total  noise  is  j~1Nj/2  +  No/2 
when  partial-band  interference  is  present  and  No/2  otherwise.  If  the  equivalent  noise 
bandwidth  of  each  bandpass  filter  in  both  the  conventional  and  the  self-normalized 
FFH/MFSK  demodulators  is  B  Hz,  then  for  each  hop  the  signal  is  received  with  noise 
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Figure  1.3:  Thermal  Noise  and  Partial-Band  Jamming  Noise  Model 


of  power  N0B  with  probability  1  —7  when  partial-band  interference  is  not  present  and 
with  noise  of  power  (7-1  Nj  +  N0)B  with  probability  7  when  partial-band  interference 
is  present. 

The  symbol  rate  is  Rs  =  Rb/  log2  M  where  M  is  the  order  of  the  MFSK  modula- 
tion and  Rb  is  the  bit  rate.  Since  the  MFSK  signal  is  assumed  to  perform  L  hops  per 
symbol,  the  hop  rate  is  Rh  =  LRS.  The  equivalent  noise  bandwidth  of  the  bandpass 
filters  in  each  of  the  M  channels  of  both  the  self-normalized  FFH/MFSK  receiver  and 
the  conventional  FFH/MFSK  receiver  must  be  at  least  as  wide  as  the  hop  rate,  and 
in  this  thesis  we  use  B  =  Rh-  The  overall  system  bandwidth  is  assumed  to  be  very 
large  compared  to  the  hop  rate. 

D.      Fast  Frequency-Hopping 

There  are  two  types  of  frequency-hopping:  fast  frequency-hopping  and  slow- 
frequency  hopping.  In  fast-frequency  hopping  the  carrier  frequency  changes  (hops) 
more  than  once  per  data  symbol  and  all  hops  within  a  symbol  duration  are  combined 
in  the  receiver  to  provide  the  sample  statistic,  but  in  slow  frequency-hopping  the 
carrier  frequency  changes  only  once  per  data  symbol.  Hence,  the  terms  of  'slow' 
and  'fast'  do  not  mean  absolutely  how  fast  the  carrier  frequency  changes.  We  define 
diversity  of  order  L  as  the  case  where  each  Af-ary  symbol  is  transmitted  L  times  over 
the  channel.  Thus,  if  Eh  is  the  energy  for  each  trasmission  then  each  M-ary  symbol 
requires  a  total  energy  of 

Es  =  LEh  (1.3) 

where  Es  is  the  symbol  energy.  The  diversity  factor  is  a  key  point  of  fast  frequency- 
hopping  communication  systems,  and  the  object  of  this  thesis  is  to  determine  under 
what  conditons  fast  frequency-hopping  systems  can  provide  immunity  to  partial-band 
interference. 


II.  ANALYSIS 

A.      Determination  of  the  Probability  of  Bit  Error 
1.      Conventional  Receiver 

An  analysis  of  the  conventional  FFH/MFSK  receiver  requires,  to  begin 
with,  a  description  of  the  samples  at  the  output  of  the  M  quadratic  detectors  be- 
fore diversity  combining.  The  partial-band  interference,  when  present,  is  assumed 
to  be  present  in  all  channels.  Assuming  equally  likely  M-ary  symbols,  we  have  the 
probability  of  symbol  error  given  L  independent  hops  as 


P°  =  %{Li  )7t'(l-7)L-t>,W  (2-1) 


where  Ps(i)  is  the  conditional  symbol  error  probability  given  that  i  out  of  L  hops 
are  jammed.  Assuming  the  signal  is  present  in  branch  one,  we  have  the  probability 
density  function  of  the  random  variable  Xik  representing  the  output  of  the  quadratic 
detector  [Ref.  3]. 

fxlk(xlk\ak)  =  ^exp  ^ j-5— J  I0  \——)u(Xih)  (2.2) 

where  ylak  is  the  signal  amplitude,  a\  is  noise  power  in  hop  k,  and  Io{*)  is  the 
zeroth-order  modified  Bessel  function  of  the  first  kind.  The  unconditional  probability 
density  function  is  found  by  evaluating 

TOO 

fxlk(xik)  =  /      fxlk{xik\ak)fAk{ak)dak  (2.3) 

The  amplitude  of  the  signal  is  modeled  as  a  Rician  distributed  random  variable  with 

a  probability  density  function  as  given  by  (1.1).    Substituting  (1.1)  and  (2.2)  into 

(2.3),  we  have 

t      t       \  1  (      xlk-\-2a2k   \       (aky/2xlk\ 

fx»{Xlk]  =  2(^  +  2^)  6XP  {-2(*l  +  2*>))  l0  [*lTwpXlk)  {2A) 
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The  probability  density  function  of  each  of  the  M  —  1  channels  with  only  Gaussian 
noise  is 

1  /        T_,.\ 

(2.5) 


•£  mk 


which  is  easily  obtained  by  replacing  the  signal  amplitude  a^  with  0  and  replacing  X\k 
with  xmki  where  m  =  2,3,4....M.  Next,  the  L  samples  are  combined  to  obtain  the 
final  sample  statistic.  Assuming  i  out  of  L  hops  are  jammed,  we  obtain  the  decision 
variable  for  branch  1 


A'i    =    £* 


A:=l 


Ik 


k=\  k=i+l 


(2.6) 


where  the  superscript  '^  implies  a  jammed  hop  and  ^  implies  a  non-jammed  hop. 

Since  each  hop  is  assumed  to  be  independent,  we  have 

®(L-0 


fxdxi)  = 


(1) 


/rtoWl') 


®    /roM?) 


ID 


(2.7) 


(2) 


where  /v(n(x^.  )  is  the  probability  density  function  when  jamming  is  present,  fy{^){x\k) 
when  no  jamming  is  present,  ®  implies  convolution,  ®i  implies  an  i— fold  convolu- 
tion, and  ®{L  —  i)  implies  an  (L  —  z)-fold  convolution.  Similarly,  for  the  noise  only 
branches,  we  have  the  random  variables 


■•'■m  /  J  A-mk 


k=l 


(2) 


k=l  k=i+l 


m  =  2, 3,  •  •  • ,  M 


and  the  probability  density  function 


/vd)(a- 


(i) 


'  yy''  V-mfc. 

TTl  k 


(2) 


/*(»>(<*) 


fxm(xm)  = 

For  MFSK,  the  probability  of  symbol  error  is  [Ref.  5] 

^(0  =  1-    r  fxi(Xi\i)\  r  fXm(Xm)dx 

JO  L-Mi 


-\®(L-t) 


M-l 


f/.;-. 


(2.8) 


(2.9) 


(2.10) 
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Equation  (2.10)  requires  fx1  {x\\i)  and  fxm(xm)-  In  order  to  obtain  /xi(a;i|i), 
we  first  rewrite  (2.7) 

/jfi(*il*)  =  /i1(*ilO®/x1(*ilO  (2-11) 

Two  probability  density  functions  are  computed  separately 

&(«.io  =  [fM\*f 

=    -C_1{K.W]'}  (2-12) 

where  Fxlk(s)  is  Laplace  transform  of  fxlk{xik\i)  and  C~l  is  the  inverse  Laplace 
transform  operation. 


njs 


exp 


-2»Q^J-(a;j-  +  2a2) 
(1  +  25«-  +  2<r*))(«rfc  +  2<r>) 


1 


(2.13) 


l+25(<7^+2<72)_ 

where  cr2.  ■  represents  jamming  noise  power  plus  thermal  noise  power.    The  inverse 
Laplace  transform  of  (2.13)  yields 


/ijU'ilO     =     exP 


2ia2k  +  xi 


2«  +  2a2)JL         i«\ 


x1{a2k  +  2a2) 


_  i^ 


x 


L2K+2^)J 


l±i 

2 


Ii-1 


°h  + 2cr2 


(2.14) 


Similarly 


/JkO&ilO     =    exp 


2{L  -ijal  +  xx 
~2^Ff2^T 


i-i-i 


2 


2(<72  +  2<r2) 


/l- 


t-i 


(L-?K    J 


al  +  2a2 


(2.15) 


where  a\  represents  thermal  noise  power.    In  order  to  obtain  /xj(xi|z),  (2.14)  and 
(2.15)  must  be  convolved.  In  general,  this  cannot  be  done  analytically. 
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For  the  noise  branch,  the  probability  density  function  is 


°-°h 


fXm(Xm) ■ 

The  Laplace  transform  of  (2.16)  is 

FxM  = 


1      --^~ 
e       kJ 


:■,< 


<lc\ 


e       >= 


[L-x) 


(2.16) 


tfJ 


s  +  8 


(i) 


■\ 


2) 


s  +  ti 


2) 


where 


tf}  =  A",      tf*^ 


H 


2a2, 


(2.17) 


(2.1s) 


Equation  (2.17)  can  be  expanded  in  a  Heaviside  partial  fraction  expansion  to  obtain 

Ku 


Fxm(s)    = 


+ 


A"n  A 12 

A'2i 


+  •••  + 


+ 


•22 


A'2(L-t) 


(2.19) 


s + /r   (s + 42^ 

where  the  coefficients  A'ij,  j  =  1,2,  •  •  •  ,i  and  A'lfc  ,  A:  =  ;  1,2,  •  •  ,L  —  i  depend  on 
L  and  i.  Specific  coefficients  for  several  values  of  L  and  i  are  given  in  Table  2.1. 
Applying  the  inverse  Laplace  transform  to  (2.19),  we  obtain 


,d) 


fxjxm)     =    e~p'<  Xm     Ku  +  K12xm  + 


2! 


+  ---  +  AW 


x 


t-i 


O'-l)!. 


?(2). 


A 


23-T 


+c-"*  x™    A'21  +  K22xm  +  -^-=M 


91 


,L-t-l 


+  A 


2(L-0 


(2.20) 


We  can  now  compute  the  second  integral  in  (2.10)  using  (2.20)  with  the 


result 


L3 


L,  i 

Coefficients 

L  =  1,    2  =  0 

#2i  =  42) 

L  =   1,      2   =   1 

#n  -  tf J 

L  =  2,    2  =  0 

#21   =  0 

#22  =  (/?[2))2 

L  =  2,    i  =  1 

A'21  =  -A'n 

L  =  2,    2  =  2 

#ii  =  o 
A'12  =  (/#>)' 

L  =  3,    2  =  0 

#21  =  0 
A'22  =  0 

#23  =  (tf })3 

L  =  3,    2  =  1 

K       -olDf        tf>        V 

#ii-P*   v-tfWv 

#21   =  tf  >tf )  -  #H^  -  ft 

1  =  3,    2  =  2 

^11    ~  Pk     Pk       -  ^2\Pk       -   7(2) 

Pk 

L  =  3,    2  =  3 

#11  =  o 
#12  =  o 

A'13  -  (A1  })3 

TABLE  2.1:  Coefficients  of  Heaviside  Partial  Fraction  Expansion 
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/     fxm{xm)dxm    =    ^— 77 
Jo  p=i  W 


ip_ 


?(i). 


P=i 


t'-l -r 


-     V  A',  e"^  Xl  T - 

£     P  13  (<  - 1  -  r)Kft(i))r+1 

L-i       r' 

A2, 


P=l 

L-t 


E^ 


?(2), 


2p 


e  ^<= 


PL"x 


'E 


.L-t'-l-r 


3=1 


^(X-i-l-rJK^V 

So  the  final  form  of  Ps{i)  which  must  be  evaluated  numerically  is 

2iot\  +  xl 


Ps(i)    =    1 


f 
Jo 


exp 


2«  +  2< 


x1(a^  +  2(72)\- 


^1- 


2 


/.-I 


exp 


2«  +  2a2) 

2(L-i)al  +  x1\  (xl{a2k 

2K2  +  2^)     , 


<  +  2<r2 


L-i-1 


1 


2 


■ 


2{al  +  2a>) 
A'iP 


Il-i-\ 


(L-i)al    J 
y/2(L  -  iWkxx 


al  +  2cr2 


E(4>' 


?n). 


P=i 


.  1  —  1  —  7" 


^  "  E  r 


p=l 


+  £ 


2p 


(2) 


=  1  (ft 
p=l 

*E 


P=i 
L-i 

P=i 


=0(r-i -ry.i^y^ 


-#>*! 


L-i-l-r 


tM-1 


dxi 


(2.21) 


(2.22) 


^(I-i-l-WiV1, 
2.      Self-normalization  Receiver 

The  difference  between  the  self-normalization  receiver  and  the  conventional 
receiver  is  the  self-normalization  factor  which  is  the  sum  of  all  M  channel  outputs 
after  sampling  at  the  outlet  of  the  quadratic  detectors.  The  new  decision  statistic  is 
the  combination  of  self-normalized  quadratic  detector  outputs.    In  order  to  find  the 
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probability  density  function  of  the  random  variables  representing  the  self-normalized 
quadratic  detector  outputs,  we  first  define  the  random  variable 

Vk  =  (X2k  +  X3k  +  •  •  •  +  XMk)  -  Xmk  ,        (m  ^  1)  (2.23) 

The  probability  density  function  of  the  sum  of  independent  random  variables  is  the 
convolution  of  the  probability  density  functions  of  those  random  variables.  Thus,  the 
probability  density  function  of  the  random  variable  Vk  requires  M  —  2  convolutions: 

fvk(vk)    =    fx2k(x2k)  ®  fx3k(x3k)  ®  •  •  •  ®  fx,k{xik)  ® 

■■■®fxMk{xMk)     (i/m)  (2.24) 

Assuming  that  all  the  probability  density  functions  of  the  branches  containing  only 
noise  are  the  same,  we  have 

fvk(vk)  =  lfxJxlk)f{M-2)       {l^l    or   m)  (225) 

The  Laplace  transform  of  fx,k{xik)  is 

^>(S)=2^7TT  (2'26) 

Hence 

/         -.  \  M-2 

2(7^7+1 

Taking  the  inverse  Laplace  transform  of  (2.27),  we  get 


*v»(3)=(0_,   ;  i  1  (2-27) 


As  can  be  seen  from  Fig.  1.2,  the  random  variables  that  represent  the  self- 
normalized  outputs  of  each  of  the  M  channels  after  diversity  combining  are  obtained 

as 

L 

£/t  =  ££W  =  l,2,...,M  (2.29) 

fc=i 
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where  Uik,  i  —  1,2,  ...,M,  h  =  1,2, . . . ,  L  are  the  random  variables  that  repre- 
sent the  self-normalized  outputs  of  each  of  the  M  channels  for  hop  A:  of  a  symbol. 
These  random  variables  are  expressed  in  terms  of  the  independent  random  variables 
Xik  ,i  =  1,2,...,  M,  k  =  1, 2, . . . ,  L  representing  the  M  independent  outputs  of  the 
M  quadratic  detectors  for  each  channel  of  the  receiver  for  hop  k  of  a  symbol  as 

Uik  =  x    +  x  Xik     +x  (2-30) 

Ait  +  A2fc  +  •  •  •  +  AmAt 
where  0  <  Uik  <  1  since  0  <  Xik  <  oo ,  z  =  1,2,...,  M.   Using  2.29,  we  can  express 
Pr(^  <  Um\i)  as 

Pr  (Ui  <  Um\i)  =  Pr  fa  Ulk  <  J2  Umk\i)  =  Pr  fa  Ulk  -  Umk  <  0\i)         (2.31) 

\k=l  k=\  /  \k-\  J 

Defining  the  random  variable  Z\k  —  U\k  —  Umk,  we  obtain 


Pr(tf!  <Um\i)    =    Pr(YZik<0\i) 


=    Pr{Z1  <0|t)  (2.32) 

where 

L 

Zl  =  £  Z1Jt  (2.33) 

fc=i 

is  our  alternative  decision  variable. 

From  2.30  and  the  definition  of  Z\k  above,  we  have 

Zk     =     U\k  —  Umk 

Xik  —  Amjt 


Aiifc  +  X2k  +  •  •  •  +  X\fk 
Xik  —  Xmk 


(2.34) 


X\k  +  Xmk  +  K- 

where  —  1  <  Zk  <  1.  This  random  variable  represents  the  difference  between  the 
power  in  the  branch  containing  the  signal  and  the  power  of  one  noise  only  branch 
which  are  detected  by  quadratic  detectors  and  normalized  by  the  self-normalization 
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factor.  We  use  the  auxiliary  variable  method  [Ref.  7]  to  obtain  the  probability  density 
function  of  the  new  random  variable  Zu-,  where  we  define  the  two  auxiliary  variables 

Wk  =  Vfc,  n  =  Xtk  (2.35) 

Reforming  (2.34)  by  replacing  Vk  and  A'n-  with  (2.35),  we  have 

Y      n(i  -  zlk)  -  zlkwk 

Xmk  = ;    „ (2.6b) 

1  +  6\k 

and  the  Jacobian  of  the  transformation  is 

J=a±^l  (237) 

2yk  +  wk 

Now  the  joint  probabilit}'  density  function  for  the  random  variables  Z\k,  Wk,  and 
Yk  is  obtained  in  terms  of  the  joint  probability  density  function  of  the  independent 
random  variables  X\k,  Xmk,  and  Vk  as 

fzlkYkwk{zlk,yk,wk\i)     =     Vj:fxlk{yk\i)fxmk\  j~— J  fvk{wk) 

2yk  +  t^fc 
(1  +  ^ife)2 

1  /       yk  +  2a\    \       (  aky/2y~k  \ 

X2(a|  +  2^)eXp^    2(*l  +  2**))Io{*l  +  2**) 
1  /    yfc(l  -  gifc)  +  znwk\ 

exP  \-in\  (2-38) 


,2^7  (M-3)!  -V  2a2. 
Before  evaluating  the  probability  density  function  of  the  random  variable  zik,  we  must 
be  careful  in  handling  the  range  of  the  variables.  Note  that  the  range  of  the  variable 
Xmk  is  —  oo  <  Xmk  <  oo,  but  recall  that  fxmk(xmk)  is  nonzero  only  for  the  range 
Xmk  >  0.  For  the  range  of  —1  <  Z\k  <  0,  Xmk  is  always  greater  than  zero;  but 
for  the  range  of  0  <  Z\k  <  1,  the  range  of  Xmk  is  —  oo  <  Xmk  <  oo.   Consequently, 
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special  care  must  be  taken  in  evaluating  fzlk{zik)-  The  probability  density  function 
of  the  random  variable  Z\k  is  obtained  from 

fOO        TOO 

fzlk(zik\i)=  /      YjrfzlkYkwk{zik,ykiWk\i)dwkdyk  (2.39) 

Jo     Jo      \J  \ 

when  —  1  <  Z\k  <  0.  Since  fxmk(*)  is  nonzero  only  for  positive  values  of  its  argument, 
then 


Jo     Jo 


oo     /,3/fc(l-Zik)M*      1 


1^1 


fzlkYkwk{zik,yk,Wk)dwkdyk  (2.40] 


when  0    <    Zu-    <    1.     Evaluating  the  integrals  in  (2.39)  and  (2.40),   we  have  the 
probability  density  function  of  random  variable  Z\k  conditional  on  i 


f°»M)  1  2  +  6d-,u)     ^  U  +  &(l-zlt) 


M-2  + 


2(1+6) 


2  +  &(l-*i*) 

-^(«ifc)u(«iib) 


1  + 


P*(l  +  *!*) 


(1 +6(1 -*i*))(l +W 


(2.41) 


where 


S^ifc)     = 


Zik(l-zlk)M-2(l+tk) 


\M-2 


exp 


-/>fc(l  ~  Z\k) 
(1  +  6(1  -   ~u-))M"1(l  +  zlk)  ~"r  Vl  +  6(1  -  *u) 

C  1  /  \  M-2 

/  PfcCifc  \ 

X 


v(M-3)!  V(l-6(l-^u))(l  +  6) 
A£^3  (   (M  -  2)(M  -  2)! 


^0  l(M-2-m)!(m!)2 
Pfc£* 

(i  +  6(i-~u-))(i+6) 

2zlk 

(1  -  zlk)(M  -  3  -  m)\ 


M-3-m  M  _  3  _  m 

+  E 


+ 


1  +  =lk 


1+6(1-^: 


M-2-m 


^0      (M  -  3  -  777  -  q)\(q\)2  V(l  +  6(1  -  *i*))(l  +  6) 
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X 


+ 


2zlk{M  -  2  -  m)2 
(l-zlk)(M  -2-m-q) 
(1  +  &(1 -*i*))(l +  *ifc)(M- 2) 

(i-*ifc)(i  +  W 


(2.42) 


and 


P* 


a! 


0" 


2a2 
6      _       °l 


(2.43) 


are  the  signal-to-noise  ratios  of  the  direct  signal  component  and  the  diffuse  (or  faded) 
signal  component,  respective^. 

Having  obtained  the  probability  density  function  of  the  random  variable 
Z\ki  we  need  to  consider  diversity  to  get  the  overall  output  sample  statistics.  The 
random  variable  Z\  is  defined 


Zi  =  X> 


u- 


fc=i 


fc=l  fc=!  +  l 

and  the  probability  density  function  of  the  random  variable  Z\  is  obtained  as 

®(L-t) 


(2.44) 


fzM\i) 


f*>MH)       ®   /*»(*!?) 


(2.45) 


It  is  difficult  to  evaluate  equation  (2.45),  so  we  leave  this  problem  for  numerical 
analysis. 

Now  the  probability  of  symbol  error  of  the  self-normalized  FFH/MFSK 


receiver  is 


P.(i)    =   Pr(U,  <  U2  1J  lh  <  U3  |J  •  •  •  U  Ui  <  UM\i) 


(2.46) 


For  M  >  2,  an  analytic  solution  of  (2.46)  is  not  possible  because  the  random  variables 
U\,U2,  Us,  •  •  •  Um-\  are  not  independent.   Hence,  we  use  the  union  bound  to  obtain 
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an  upper  bound  on  system  performance.  From  (2.46) 

Ps(i)    <    Pr(Ui  <  U2\i)   +    Pr(Ux  <  U3\i)   +  ••■  +   Pr(Ux  <  UM\i)  (2.47) 

Since  we  assume  the  sample  statistics  of  all  noise  only  branches  are  the  same,  then 
(2.47)  becomes 

Ps(i)    <    (M  -  1)  Pr{Ux  <  Um\i)  (2.48) 

where  m  =  2, 3, 4,  •  •  • ,  M.  Substituting  (2.32)  into  (2.48),  we  have 

Ps(i)  <  (M  -  l)Pr{Zx  <  0|0  (2.49) 


Since 


then  from  (2.49) 


Pr{Zx  <0)    =    /     fzAzl\i)dzx  (2.50) 

Ps(i)   <   (M-l)  f°  /Zl(*i|i)«fei  (2.51) 

The  probability  of  symbol  error  is  related  to  the  probability  of  bit  error  by  [Ref.  5] 

cyk-\ 

Pb(i)     =     Ps(i)    x    2F3Y'*   =   l°9^M 

=  ™  x  m^T)  (2-52) 

Substituting  (2.48)  into  (2.52),  we  get 

Pb(i)    <    y  fLfzM\i)dzx  (2.53) 

Equation  (2.53)  must  be  evaluated  numerically. 
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III.  NUMERICAL  ANALYSIS 

Results  of  the  probability  of  error  analysis  for  both  the  conventional  receiver 
and  the  self-nomalization  receiver  are  obtained  for  L  =  1  to  6  hops  per  bit  for  several 
signal-to-noise  ratios  and  direct-to-diffuse  component  ratios  as  a  function  of  signal-to- 
interference  ratio  to  provide  an  overall  comparison  of  channel  and  interference  effects 
on  system  performance  for  both  receivers.  In  all  cases  the  evaluation  of  either  (2.22) 
or  (2.41)  is  accomplished  via  a  numerical  integration  routine.  We  consider  L  +  1 
different  cases  for  a  diversity  of  L  hops.  For  example,  for  L  =  1  the  two  cases  to  be 
considered  are  one  non-jammed  hop  and  one  jammed  hop.  For  L  =  2  there  are  three 
cases:  neither  hop  jammed,  one  non-jammed  hop  and  one  jammed  hop,  and  both  hops 
jammed.  Additionally,  the  effect  of  the  jamming  ratio  7  on  system  performance  is 
examined.  We  can  relate  the  signal-to-noise  ratio  to  the  signal  energy  density-to-noise 
power  ratio  for  a  jammed  hop  as 

a2  +  2a2  Eb         B  log2  M 


where 


<  No  +  NjhR,      L 

a\  :  direct  signal  power  component 

2a2  :  diffuse  signal  power  component 

a2,  -f  2<j2    :  total  signal  power 

a2.^  :  thermal  noise  plus  jamming  noise  power 

Eb  :  average  signal  bit  energy  density 

./Vo  :  thermal  noise  power  spectral  density 

Nj  :  jamming  noise  power  spectral  density 

7  :  jamming  ratio 
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(3.1) 


B  :  IF  bandwidth 

Rs  :  symbol  rate 

M  :  number  of  symbols  in  alphabet 

L  :  number  of  hops  in  a  symbol 

For  hops  with  no  partial-band  jamming 

a\  +  2cr2  Eb_  B  log2  M 

<j\         "    N~0R~S      L 

where  a\  is  the  thermal  noise  power.  Assuming  the  IF  bandwidth  (B)  is  equal  to  the 

symbol  rate,  we  can  rewrite  (3.1)  and  (3.2) 

cc\  +  2a2  Eh         log2  M 


(3.2) 


and 

a\  +  2cr2  Eb  log2  M 

o\  N~0      L 

For  computational  purposed,  define 


akj 

— 

1 

Eb 

No 

— 

V 

Eb 

Nj 

V 

2a2 

= 

i> 

log2  M 

= 

W 

Substituting  (3.5)  through  (3.9)  into  (3.3)  and  (3.4),  we  have 


(3.3) 
(3.4) 

(3.5) 
(3.6) 

(3.7) 

(3.8) 
(3.9) 


2  rnvij)  u 

k  O7  +  7")(0  +  l)£ 

9  7^ 

°l      =      7-T—T  3'12 

(77  +  7^) 
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Substituting  (3.10),  (3.11),  and  (3.12)  into  the  probability  density  function  of  the 
conventional  receiver  and  the  probability  density  function  of  self-normalization  re- 
ceiver and  using  an  fast  Fourier  transform  (FFT)  routine  with  a  sequence  length  of 
N=1024,  numerical  sequences  of  the  probability  density  functions  are  generated  for 
various  signal-to-noise  ratios  of  the  jammed  and  non-jammed  cases. 


A.      Conventional  Receiver 

Except  for  a  single  convolution,  we  have  obtained  analytic  solution  for  the  prob- 
ability density  function  of  the  decision  statistic  for  the  conventional  FFH/MFSK  de- 
tector. Unfortunately,  the  modified  Bessel  function  terms  make  it  difficult  to  obtain 
a  complete  analytic  solution  for  the  general  case.  Consequently,  we  use  an  FFT  to 
perform  the  convolution  numerically.  Even  though  we  have  an  analytic  expression  for 
the  probability  density  function  of  the  noise  branch  statistic  (the  last  part  of  (2.22)), 
we  use  numerical  methods  to  correspond  with  the  signal  statistic  part  which  must  be 
evaluated  numerically.  Hence,  we  use  different  equations  for  the  numerical  procedure. 
Equation  (2.16)  can  be  rewritten 

fxm(xm)  =  f[xl{Xm)    ®    f[xl(Xm) 


where 

JXm\Xm)        = 
fxm{Xm)        = 

Using  the  Laplace  transform,  we  get 


*&« 


e       ki 


2°h 


1  -*± 

e       * 

2*1 


®(L-t) 


(3.13) 

(3.14) 
(3.15) 


ft> 


s  +  Pl1] 


0 


(2) 


-i  L-i 


s  +  tf ) 


(3.16) 
(3.17) 
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where 


tfJ 


i(2>         - 


?VT? 


la 


fcj 


/9l  ;      -      — — 


2*2 


CMS) 
(3.19) 


The  inverse  Laplace  transform  of  (3.16)  and  (3.17)  yields 


fx,n(xm) 


fxm(Xm)      ~ 


1 


2< 
23. 


.«'-! 


(t-1) 


fcj 


L-i 


,L-i-\ 


■e    *% 


(3.20) 
(3-21) 


>k/      (L-i-iy. 

Now  we  have  a  different  form  of  the  probability  density  function  of  the  random  vari- 
able Xm 

i  V 


fxm{Xr 


2ah        (t-1)! 


t-i       _^m_ 
e      *■* 


'fcj 


•e      * 


(3.22) 


(L-t-1)! 

Instead  of  (2.21)  we  substitute  (3.22)  into  (2.22)  and  analyze  (2.22)  numerically  be- 
cause we  need  a  numerical  sequence  of  this  formula  to  coincide  with  the  numerical 
sequence  obtained  for  (2.11)  The  length  of  the  basic  sequence  is  N=1024,  and  prior 
to  convolution  it  is  zero-padded  to  create  a  sequence  of  total  length  N=2048  since  the 
sequence  is  to  be  convolved  one  time. 

Numerical  problems  are  encountered  since  as  the  value  of  the  variable  X\  in- 
creases the  exponential  terms  approach  zero  and  the  modifiied  Bessel  function  terms 
approach  infinity.  Hence,  the  multiplication  of  the  exponential  term  and  the  modified 
Bessel  function  term  is  invalid  when  the  argument  of  the  exponent  is  less  than  -745 
or  the  argument  of  the  modified  Bessel  function  is  greater  than  705  when  386Mat- 
lab  is  used.  For  example,  in  the  situation  described  above,  the  value  of  exp(-745)  is 
4.9470e-324,  but  the  value  of  exp(-746)  is  set  to  0  by  the  computer.    The  value  of 
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Bessel(705),where  the  386  MATLAB  function  for  the  modified  Bessel  function  Io{*) 
is  Bessel(»),  is  2.2621e+304  and  Bessel(706)  is  not  defined.  Even  though  the  product 
of  the  two  is  either  zero  or  an  invalid  number,  the  true  value  of  a  multiplication  is  a 
valid  number.  When  this  situation  is  encountered,  the  program  evaluates  the  situa- 
tion and  uses  either  zero  or  one  (because  the  cumulative  probability  density  function 
approaches  one  as  the  value  of  the  variable  increases). 

When  two  probability  density  functions  are  convolved  by  an  FFT,  if  the  valid 
range  of  one  probability  density  function  is  much  wider  than  the  range  of  the  other,  the 
result  is  almost  the  same  as  the  probability  density  function  that  has  wider  range. 
We  must  check  for  this  situation  whenever  the  parameters  of  the  equation  change 
while  the  program  is  running,  skip  the  numerical  convolution  steps,  and  instead  use 
the  probability  density  function  that  has  wider  range  as  the  convolution  output. 

B.      Self-Normalization  Receiver 

The  probability  density  function  of  self-normalization  receiver  (2.41)  is  relatively 
complicated  and  we  must  use  an  L-io\d  convolution  to  evaluate  (2.44);  an  i-fold  con- 
volution for  jammed  hops,  an  (L  —  z)-fold  convolution  for  nonjammed  hops,  and  each 
result  must  be  convolved.  To  implement  the  convolution  numerically,  the  numerical 
sequences  of  the  probability  density  function  of  the  self-normalization  receiver  must 
be  zero-padded  to  a  total  length  of  N  x  L  before  the  FFT  is  performed  in  order 
to  preserve  the  desired  linear  convolution  of  (2.44).  The  transformed  sequences  are 
multiplied  point  by  point  to  produce  the  desired  seqence  of 

**w  =  [*£wf  x  \FziSs)]L~l  (3-23) 

where  F  is  the  Fourier  transform  of  /.  This  result  is  then  inverse  fast  Fourier  trans- 
formed to  yield  the  desired  sequence  for  fzi{z\).  The  procedures  for  L=3  and  L=b 
are  similar  to  that  for  L=A  and  L=6,  respectively,  because  the  length  of  the  sequence 

26 


for  the  FFT  should  be  a  power  of  two,  and  the  sequences  should  be  zero-padded  to 
7V=4096  (212)  for  L=3  and  L=4  and  to  JV=8192  (213)  for  L=5  and  L=6. 

Equation  (2.42)  has  many  summations  inside  which  must  be  carefully  handled 
since  numerical  errors  are  likely  when  the  value  of  the  variable  Z\k  approaches  one 
When  Z\k  approaches  unity,  the  term  1/(1  —  Z\k)  approaches  infinity.  Therefore,  we 
must  group  the  terms  properly  to  minimize  numerical  errors.  Equation  (2.42)  is 
rewritten 


9(zik)     = 


iM-2 


^(i-^u-r-2(i  +  6 

(i  +  6(i-^u))M-1(i  +  ^ 

^2   (M-2)(M-2)! 


exp 


M  -3 

+  Yl  exp 

771=0 


'0  {M  -2-m)!(m!)2 


-pk{l  -  Zlk) 

l  +6(1  -2ut) 

PkZlk 


(i  +  6(i-^u))(i  +  6) 


i  +  6(i-^u-). 
2,2,((1  -  zlk)(l  +  6))M"3-m(l  +  6)(1  +  zlkr~\M  -  2  -  m) 


(l+ai-^))^1-"1 


M-2-m 


M  -  2  - 


?/; 


g5o     (M -2 -m -*)!($!)* 

M-3 


PkZ\k 


+  H  exp 

771=0 

A/-3-771 


■ 


E 


~^(l-~l/c)    1    ((1 
l+6(l-2lfc) 

M  -  3  -  m 


(i+6(i-2u))(i+6). 

fc)(l  +  6))M"3-m^u(l  +  zn)m(M  -  2] 


-i 


g-     (M -3 -m- *)!(?!)' 

in  order  to  minimize  computational  errors. 


(l+^(l-^))M-2-, 


(i  +  6(i-^))(i  +  6)J 


(3.24) 
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IV.  RESULTS 

A  FFH/MFSK  receiver  design  that  demonstrates  strong  immunity  to  the  partial- 
band  interference  is  the  noise  normalization  receiver.  This  receiver  requires  a  measure- 
ment of  the  noise  power  per  hop,  but  it  is  not  easy  to  measure  the  noise  power  precisely 
in  practical  situations.  This  measurement  is  not  required  for  the  self-normalization 
receiver;  hence,  it  is  easier  to  implement  the  self-normalization  receiver.  As  will 
be  seen,  partial-band  jamming  strongly  affects  the  performance  of  the  conventional 
FFH/MFSK  receiver,  while  the  self-normalization  FFH/MFSK  receiver  performs  well 
against  partial-band  jamming.  Probability  of  bit  error  analyses  for  both  the  conven- 
tional system  and  the  self-normalization  system  assume  worst  case  partial-band  jam- 
ming. The  channel  is  modeled  for  a  moderately  strong  direct-to-diffuse  signal  ratio 
{a2k/ 2a2  =  10)  which  we  will  refer  to  as  Rician,  and  as  a  primarily  Rayleigh  faded 
channel  (al/2a2  =  0.01).  For  Rician  fading  the  worst  case  jamming  ratio  typically 
increases  as  interference  power  increases.  For  Rayleigh  fading  worst  case  jamming 
corresponds  to  broadband  jamming  (jamming  ratio  7=1)  for  all  levels  of  jamming 
power. 

A.      Performance  of  Conventional  FFH/MFSK  Receiver 

Graphs  of  the  probability  of  bit  error,  including  worst  case,  as  a  function  of 
signal-to-jamming  ratio  for  fixed  signal-to-thermal  noise  ratios  are  obtained  for  1 
through  6  hops  per  symbol.  These  results  are  obtained  by  numerically  evaluating 
(2.44)  for  values  of  the  jamming  ratio  ranging  from  1.0  to  0.001  and  retaining  worst 
case  performance  for  each  value  of  E^/Nj  .  The  results  for  a  Rician  fading  channel 
are  shown  in  Fig.  4.1  through  Fig.  4.6. 
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Figure  4.1:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  1,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.2:    The  Performance  of  the  Conventional   FFH/MFSK  none 
herent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  2,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.3:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  3,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.4:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.5:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  5,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.6:    The  Performance  of  the  Conventional   FFH/MFSK   none, 
herent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  6,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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The  signal-to-noise  ratio,  direct-to-diffuse  ratio,  and  the  modulation  order  are 
chosen  to  enable  direct  comparison  with  the  results  obtained  in  [Ref.  2].  We  see  that 
for  a  Rician  channel  that  diversity  is  not  effective  in  reducing  the  effects  of  worst  case 
partial-band  jamming.  Worst  case  performance  is  poorer  when  diversity  is  used  than 
for  no  diversity.  As  diversity  increases,  only  the  asymtotic  value  for  high  Eb/Nj  ap- 
proaches the  performance  of  the  noise-nomalization  receiver.  The  order  of  modulation 
has  no  effect  with  respect  to  providing  immunity  to  partial-band  interference.  This 
is  illustrated  in  Fig.  4.4  and  Fig.  4.7  through  Fig.  4.9.  Generally  both  diversity  and 
the  order  of  modulatin  are  effective  in  improving  system  performance  for  full-band 
jamming  (7  =  1)  as  shown  Fig.  4.10  and  Fig.  4.11. 

Generally  both  diversity  and  the  order  of  modulation  for  a  Rayleigh  fading 
channel  are  effective  in  improveing  system  performance  for  full  band  jamming  (7=1) 
as  shown  in  Fig.  4.10  and  Fig.  4.11.  However,  System  performance,  shown  in  Fig. 
4.12  through  Fig.  4.16,  is  not  improved  by  either  diversity  or  modulation  order  when 
worst  case  partial-band  jamming  is  assumed.  When  the  signal-to-thermal  noise 

ratio  increases,  the  performance  of  the  conventional  receiver  for  worst  case  partial- 
band  jamming  is  unchanged  from  that  of  lower  signal-to-thermal  noise  cases  except 
for  the  asymtotic  limit  of  negligible  Eb/Nj. 
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Figure  4.7:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  8  and  Direct-to-DifFuse 
Signal  Power  Ratio  =  10 
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Figure  4.8:     The  Performance  of  the  Conventional   FFH/MFSK   non< 
herent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  16  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.9:    The  Performance  of  the  Conventional   FFH/MFSK   non 
herent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  32  and  Direct-to-DifFuse 
Signal  Power  Ratio  =  10 
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Figure  4.10:  The  Performance  of  the  Conventional  FFH/MFSK  noncoh< 
ent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0dB,  Jamming  Ratio 
=  1,  Order  of  Modulation  =  8  and  Direct-to-DifFuse  Signal  Power  Ratio 
=  0.01 
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Figure  4.11:  The  Performance  of  the  Conventional  FFH/MFSK  noncoher- 
ent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0dB,  Hopping  Num- 
ber per  Symbol  =  4,  Jamming  Ratio  =  1  and  Direct-to-DifFuse  Signal 
Power  Ratio  =  0.01 
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Figure  4.12:    The  Performance  of  the  Conventional  FFH/MFSK  non. 
herent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  1,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  0.01 
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Figure  4.13:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  3,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  0.01 
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Figure  4.14:  The  Performance  of  the  Conventional  FFH/MFSK  non 
herent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  5,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  0.01 
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Figure  4.15:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  16  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  0.01 
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Figure  4.16:  The  Performance  of  the  Conventional  FFH/MFSK  nonco- 
herent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  32  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  0.01 
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B.      Performance  of  Self-Normalization  FFH/MFSK  Receiver 

All  parameters  used  in  the  graphs  for  the  self-normalization  receiver  are  the 
same  as  those  used  for  the  conventional  receiver  to  make  it  easier  to  compare  with 
each  other  as  well  as  with  the  noise-normalization  receiver.  As  we  see  from  Fig.  4.17 
through  Fig.  4.25,  for  a  Rician  fading  channel,  increasing  diversity  and  increasing 
modulation  order  are  both  effective  in  minimizing  the  effects  of  worst  case  partial- 
band  interference. 

Overall,  the  performance  of  the  self-normalization  receiver  is  almost  the  same 
as  that  of  the  noise-normalization  receiver.  There  exists  a  small  difference  in  the 
performance  of  the  two  receivers.  This  difference  is  expected  due  to  the  union  bound 
used  to  obtain  the  FFH/MFSK  self-normalized  receiver  performance.  The  union 
bound  has  the  most  pronounced  effect  on  system  performance  when  Eb/Nj  is  small 
and  the  modulation  order  is  high.  See,  for  example,  Fig.  4.26  where  the  probability 
of  bit  error  for  a  bit  energy- to- jamming  noise  density  ratio  of  0  dB  with  a  jamming 
ratio  1.0  is  shown  greater  than  unity.  This  is  not  valid  for  a  probability.  However, 
the  error  performance  at  low  Eb/Nj  is  not  important,  and  the  error  due  to  the  union 
bound  compared  to  the  exact  value  at  high  Eb/Nj  is  not  significant. 

The  performance  of  the  self-normalization  receiver  for  a  Rayleigh  fading  channel 
in  worst  case  partial-band  jamming  is  the  same  as  for  full  band  jamming  (7=1) 
regardless  of  the  parameters  as  shown  in  Fig.  4.27  through  Fig.  4.29.  Generally  for 
a  Rayleigh  fading  channel,  diversity  is  effective  in  improving  the  performance  of  the 
self-normalization  receiver.  This  is  shown  in  Fig.  4.30.  The  effect  of  increasing  the 
modulation  order  is  less  significant  as  shown  in  Fig.  4.26.  As  mentioned  previously, 
the  union  bound  makes  the  simulated  performance  worse  than  the  true  performance. 
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Figure  4.17:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  1,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.18:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  2,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.19:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  3,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.20:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  4  and  Direct-to-DifFuse 
Signal  Power  Ratio  =  10 
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Figure  4.21:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  5,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.22:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  6,  Order  of  Modulation  =  4  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.23:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  8  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.24:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  16  and  Direct-to-DifFuse 
Signal  Power  Ratio  =  10 
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Figure  4.25:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  13.35  dB,  Hopping 
Number  per  Symbol  =  4,  Order  of  Modulation  =  32  and  Direct-to-Diffuse 
Signal  Power  Ratio  =  10 
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Figure  4.26:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0dB,  Hopping 
Number  per  Symbol  =  4,  Jamming  Ratio  =  1  and  Direct-to-Diffuse  Signal 
Power  Ratio  =  0.01 
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Figure  4.27:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0  dB,  Hopping 
Number  =  1,  Order  of  Modulation  =  4  and  Direct-to-DifFuse  Signal  Power 
Ratio  =  0.01 
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Figure  4.28:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0dB,  Hopping 
Number  =  4,  Order  of  Modulation  =  4  and  Direct-to-Diffuse  Signal  Power 
Ratio  =  0.01 
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Figure  4.29:  The  Performance  of  the  Self-Normalization  FFH/MFSK  non- 
coherent Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0dB,  Hopping 
Number  =  6,  Order  of  Modulation  =  32  and  Direct-to-Diffuse  Signal  Power 
Ratio  =  0.01 
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Figure  4.30:  The  Performance  of  the  Self-Normalization  FFH/MFSK  n 
coherent  Receiver  at  Signal-to-Thermal  Noise  Ratio  =  20.0dB,  Jamming 
Ratio  =  1,  Order  of  Modulation  =  8  and  Direct-to-DifFuse  Signal  Power 
Ratio  =  0.01 
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V.  CONCLUSION 

The  probability  of  bit  error  performance  for  a  conventional  noncoherent  FFH/MFSK 
receiver  has  been  obtained  for  Rayleigh-  and  Rician-faded  channel  with  partial- 
band  jamming.  The  effect  of  fading  is  detrimental  for  both  fading  models,  but  fast 
frequency-hopping  in  general  provides  a  means  to  overcome,  at  least  partially,  fading 
effects.  However,  the  conventional  receiver  is  severely  affected  by  worst  case  partial- 
band  jamming.  It  does  not  show  any  improvement  in  worst  case  performance  with 
high  diversity  or  high  order  of  modulation.  The  probability  of  bit  error  of  the  conven- 
tional receiver  in  worst  case  partial-band  jamming  is  always  between  10~2  and  10-' 
at  Eb/Nj  of  20  dB  regardless  of  Eb/No,  the  order  of  modulation,  diversity  factor,  or 
direct-to  diffuse  power  ratio.  The  probability  of  bit  error  of  the  self-normalization  re- 
ceiver in  worst  case  partial-band  jamming  is  the  same  as  that  of  conventional  receiver 
for  slow  frequency-hopping  (L=l).  However,  for  fast  frequency-hopping,  the  perfor- 
mance becomes  better  and  better  with  increasing  diversity  factors:  the  probability  of 
bit  error  is  5  x  10"4  at  L=2,  2  x  10~4  at  L=3,  10-4  at  L=4  and  8  x  lO-5  at  L=5  and 
L=6.  When  the  diversity  factor  is  6,  the  worst  case  performance  is  almost  same  as 
that  of  full  band  jamming  case  (jamming  ratio  1.0).  At  this  point  the  jammer  does 
not  need  any  information  about  the  communicator  to  optimize  the  jamming  effect  by 
managing  the  jamming  power. 

The  self-normalization  receiver  with  diversity  provides  very  good  immunity  to 
worst  case  partial-band  jamming  and  has  a  performance  similar  to  that  of  the  noise- 
normalization  receiver.  The  order  of  modulation  also  improves  the  overall  perfor- 
mance of  the  self-normalization  receiver,  but  partial-band  jamming  affects  higher 
orders  of  modulation  more  than  the  lower  ones.    Thus,  if  the  communicator  uses  a 
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higher  order  of  modulation  and  wants  to  suppress  all  partial-band  jamming  effects, 
he  should  use  a  higher  hopping  rate. 
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APPENDIX  A 


%  This  is  the  Matlab  program  to  produce  the  probability  of  bit   % 

error  of 

%  Conventional  FH/MFSK  quadratic  receiver. 

%  Lee,Kang  Yeun   (Major  ROKAF) 

clear 
%  3  signal-to-noise  ratios  to  be  implemented  (13.35  db,  16  db  and 
%  20  db) 

SNR=[21.62718524D0  39 . 81071706D0  100. 0D0]  ; 
%  2  direct-to-diffuse  signal  power  ratios  (0.01  and  10) 

DDR=[0.01D0  10.0D0]; 
%  4  jamming  ratio  (1,  0.1,  0.01  and  0.001) 

JRATIO=[1.0D0  0.1D0  0.01D0  0.001D0]; 

data=zeros (8, 51)  ; 
%  set  the  thermal  plus  jamming  noise  power  1 

rhokl=l; 
%  basic  sequence  length  is  1024 

lim=1024; 

n=2*lim; 

zpad=zeros(l, lim) ; 

opad=ones (1, lim) ; 
%  hopping  number  (1  to  6) 

for  L=l:6, 

for  jl=l:2,   psi=DDR(jl); 

for  j 2=1:3,   nu=SNR(j2) ; 

for  j3=l:4,   gama=JRATIO(j3) ; 
%  the  order  of  modulation  (M-ary  :  4,  8,  16,  32-ary) 

for  j4=2:5;   M=2Aj4; 

for  j5=0:2:50,  iota=10. 0D0A (0. 1D0*J5) ; 
%  diffuse  signal  power  (  rho  ) 

rho=nu*gama*iota/ (nu+gama*iota)/ (psi+1) *j4/L; 
%  direct  signal  power  (  alpha  ) 

alpha=rho*psi ; 
%  thermal  noise  power  (  rhok2  ) 

rhok2=iota*gama/ (iota*gama+nu) ; 
%  set  the  sequence  limit  roughly  based  on  the  singal-to-noise  ratio 
%  (nu) 

limx=8  0*nu  ; 
%  divide  the  range  by  100  to  check  the  proper  limit 

delta=limx/100 ; 
%  select  almost  the  middle  value  of  the  increment  to  compute  the 
%  sequence 
%  of  the  probability  density  function. 
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x=0 . 49*delta : delta : limx-0 . 51*delta ; 

g=exp(-(x+2*L*alpha)/2/ (rhokl+rho) )/2/ (rhokl+rho) .  ... 
*  x/2/L/alpha) 
( (L-l)/2) .*abs(jA (L-l) *bessel (L-l, j*sqrt (2*L*alpha*x)  . . . 
/ (rhokl+rho) ) ) *delta ; 
%  check  the  sequence  whether  it  has  invalid  sequence  or  not 

fi=find(isnan(g) ) ; 

g(fi)=zeros(l, length (fi) ) ; 
%  if  the  sequence  has  only  Os,  the  limit  is  too  wide.  Thus  set  the 
%  first 
%  increment  as  the  new  limit 

if    all(g==0)==l, 

gmin=l; 
%  if  the  sequence  is  valid,  check  the  point  around  1.0e-6  to  select 
%  the 
%  new  limit 

else,  [xx,gmax]=max(g) ; 

g(l:gmax) =ones (l,gmax) *xx; 

[xx,gmin]=min( (g-1.0e-6) .  A2)  ; 

end 

1 imx=gmin*delta ; 
%  if  the  new  limit  is  first  or  second  value  of  the  increment 
%  redo  above  procedure 

while  gmin  <=  2  , 

delta=limx/100;  x=0. 5*delta: delta: limx-O. 5*delta ; 

g=exp (-(x+2*L*alpha)/2/ (rhokl+rho) )/2/  . . . 
(rhokl+rho) .* (x/2/L/alpha) . A  ... 

( (L-l)/2) . *abs(jA (L-l) *bessel (L-l, j*sqrt (2*L*alpha*x)  . . . 
/ (rhokl+rho) ) ) *delta ; 

fi=f ind(isnan(g) ) ; 

g (fi)=zeros(l, length (fi) )  ; 

if    all(g==0)==l, 

gmin=l; 

else,  [xx,gmax]=max(g) ; 

g ( 1 : gmax) =ones ( 1 , gmax) *xx ; 

[xx,gmin]=min( (g-1. Oe-6) .  A2)  ; 

end 

1 imx=gmin*delta ; 

end 
%  after  setting  the  limit,  whole  range  is  divided  by  1024 
%  to  produce  the  sequence 

delta=limx/lim;x=0.4  9*delta:delta:limx-0.51*delta; 
%  if  jamming  ratio  is  1,  the  computation  is  simpler  than  others 

if  j3  ==  1, 

gs=-(x+2*L*alpha)/2/ (rhokl+rho) -log (2* (rhokl+rho) ) 

+log(x/2/L/alpha)*  ... 

(L-l)/2+log(abs(jA (L-l) *bessel (L-l, j*sqrt (2*L*alpha*x)  . . . 

/(rhokl+rho)))) ; 

gs=exp (gs) *delta ; 

f i=f ind(isnan(gs)  )  ; 
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gs(f i)=zeros (1, length (fi) )  ; 

if  all(gs==0)==l, 

gs=zeros (1, lim) ; 

gs(l)=l; 

else,gs=gs/sum(gs) ; 

end 

ps=log (l/2/rhokl)*L+log(x)*(L-l) -log ( fact (L-l) )  -  x/2/rhokl; 

ps=exp(ps) *delta; 

fi=f ind (isnan (ps)  )  ; 

ps(fi)=zeros (1, length (fi) ) ; 

if  all(ps==0)==l, 

ps=ones(l, lim) ; 

else,ps=ps/sum(ps) ; 

ps=cumsum(ps) ; 

ps=ps/max(ps) ; 

end 

data (count, j5+l)=M/2/ (M-l) * (l-sum(gs. *ps. A (M-l) ) ) ; 

else 
%  if  jamming  ratio  is  not  1  ,  the  hopping  number  has  to  be 
%   considered 

for  i  =  0  :  L 
%  if  no  hop  is  jammed 

if  i  ==  0, 

delta=limx/100; 

x=0.5*delta: delta: limx-0. 5*delta; 

g=-(x+2*L*alpha)/2/(rhok2+rho)-log(2*(rhok2+rho) )  ... 

+log(x/2/L/alpha) *  ... 

(L-l)/2+log(abs(jA(L-l)*bessel(L-l, j*sqrt(2*L*alpha*x)  . . . 

/(rhok2+rho) )) ) ; 

g=exp(g) *delta; 

f i=f ind (isnan (g) ) ; 

g(f i)=zeros(l, length (fi) ) ; 

if  all(g==0)==l, 

gmin=l; 

else, [xx,gmax]=max(g) ;g(l:gmax)=ones (l,gmax) *xx; 

[xx,gmin]=min( (g-1.0e-6) . A2) ; 

end 

vlimx=gmin*delta; 
%  check  the  sequence  limit  again 

while  gmin  <=  2, 

delta=vlimx/100 ; 

x=0. 5*delta:delta:vlimx-0.5*delta; 

g=-(x+2*L*alpha)/2/(rhok2+rho)-log(2*(rhok2+rho) )  ... 

+log(x/2/L/alpha) *  ... 

(L-l)/2+log(abs(jA (L-l) *bessel (L-l, j*sqrt (2*L*alpha*x)  . . . 

/(rhok2+rho) ) ) ) ; 

g=exp (g) *delta ; 

f i=f ind ( isnan (g) ) ; 

g(f i)=zeros(l, length (fi) ) ; 
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if  all(g==0)==l, 

gmin=l; 

else, 

[xx,gmax]=max(g) ;g(l:gmax)=ones (l,gmax) *xx; 

[xx,gmin]=min( (g-1.0e-6) . A2) ; 

end 

vlimx=gmin*delta; 

end 

delta=vlimx/lim; 

x=0.49*delta:delta:vlimx-0.51*delta; 

gs=-(x+2*L*alpha)/2/(rhok2+rho)-log(2*(rhok2+rho) )  ... 

+log (x/2/L/alpha) *  . . . 
(L-l)/2+log(abs(jA(L-l)*bessel(L-l,j*sqrt(2*L*alpha*x)  . . . 
/(rhok2+rho)))) ; 
gs=exp (gs) *delta ; 
fi=find(isnan(gs) ) ; 
gs(f i)=zeros( 1, length (f i) ) ; 
if  all(gs==0)==l, 
gs=zeros(l,lim) ; 
gs(l)=l; 
else, 

gs=gs/sum(gs) ; 
end 

ps=-log(2*rhok2)*L+log(x)*(L-l)-log(fact(L-l) )  -  x/2/rhok2 ; 
ps=exp(ps) *delta; 
f i=f ind(isnan(ps) ) ; 
ps(fi)=zeros(l, length (fi) ) ; 
if  all(ps==0)==l/ 
ps=ones(l,lim) ; 
else,  ps=ps/sum(ps) ; 
ps=ps/sum(ps) ; 
ps=cumsum(ps) ; 
ps=ps/max(ps) ; 
end 

if  all  hops  are  jammed 
elseif  i  ==  L, 
delta=l imx/ 100 ; 

x=0 . 5*delta : delta : limx-0 . 5*delta ; 

g=-(x+2*i*alpha)/2/(rhokl+rho)-log(2*(rhokl+rho) )  — 
+log(x/2/i/alpha)*  ... 

(i-l)/2+log(abs(jA(i-l)*bessel(i-l, j*sqrt(2*i*alpha*x)  . . . 
/(rhokl+rho)))) ; 
g=exp(g) *delta; 
f i=find(isnan(g) ) ; 
g (fi)=zeros(l, length (fi) ) ; 
if  all(g==0)==l, 
gmin=l ; 
else, 

[xx,gmax]=max(g) ;g(l:gmax)=ones (l,gmax) *xx; 
[xx,gmin]=min( (g-1.0e-6) . A2) ; 
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end 

ulimx=gmin*delta; 
check  the  sequence  limit  again 
while  groin  <=2, 
delta=ulimx/100; 

x=0.5*delta: delta :ulimx-0.5*delta; 

g=-(x+2*i*alpha)/2/ (rhokl+rho) -log (2* (rhokl+rho) )  ... 
+log(x/2/i/alpha) *  ... 

(i-l)/2+log(abs(jA(i-l)*bessel(i-l, j*sqrt (2*i*alpha*x)  . . . 
/(rhokl+rho))) ) ; 
g=exp (g) *delta ; 
f i=f ind (isnan(g) ) ; 
g(f i)=zeros (1, length (fi) )  ; 
if  all(g==0)==l, 
gmin=l; 
else, 

[xx/gmax]=max(g) ; 
g(l:gmax)=ones (l,gmax) *xx; 
[xx, groin ]=roin( (g-1. Oe-6) . ^2) ; 
end 

ulimx=groin*delta ; 

end 

delta=ul imx/1 im ; 

x=0 . 49*delta : delta : ulimx-0 . 51*delta ; 

gs=-(x+2*i*alpha)/2/ (rhokl+rho) -log (2* (rhokl+rho) )  ... 
+log(x/2/i/alpha) *  ... 
(i-l)/2+log(abs(jA(i-l)*bessel(i-l, j*sqrt (2*i*alpha*x)  . . . 

/(rhokl+rho)) )) ; 

gs=exp(gs) *delta; 

f i=f ind(isnan(gs) ) ; 

gs(f i)=zeros (1, length (fi) ) ; 

if  all(gs==0)==l, 

gs=zeros(l, lim) ; 

gs(l)=l; 

else, 

gs=gs/sum(gs)  ; 

end 

ps=-log(2*rhokl)*i+log(x)*(i-l)-log(fact(i-l) )  -  x/2/rhokl; 

ps=exp (ps) *delta ; 

if    all(ps==0)  ==1, 

ps=  ones(l, lim) ; 

else, 

fi=  f ind (isnan(ps) ) ; 

ps(f i)=zeros(l, length (fi) ) ; 

ps  =  cumsum(ps) ; 

ps  =  ps/max(ps); 

end 
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else 
%  if  i  hops  are  jammed  and  L-i  hops  are  not  jammed 
%  check  the  sequence  limit  again 

rhok=rhokl ; 

delta=limx/100; 

x=0 . 5*delta : delta : 1 imx-0 . 5*delta ; 

g=-(x+2*i*alpha)/2/ (rhok+rho) -log (2* (rhok+rho) )  ... 

+log (x/2/i/alpha) *  . . . 

(i-l)/2+log(abs(jA(i-l)*bessel(i-l, j*sqrt (2*i*alpha*x) 

/(rhok+rho)))) ; 

g=exp (g) *delta ; 

f i=find(isnan(g) ) ; 

g (fi)=zeros(l, length (fi) ) ; 

if  all(g==0)==l, 

gmin=l; 

else, 

[xx, gmax]=max(g) ; 

g ( 1 : gmax) =ones ( 1 , gmax) *xx ; 

[xx,gmin]=min( (g-1.0e-6) .  A2) ; 

end 


wlimx=gmin*delta; 

while  gmin  <=2 , 

delta=wlimx/100 ; 

x=0 . 5*delta : delta : wl imx-0 . 5*delta ; 

g=-(x+2*i*alpha)/2/ (rhok+rho) -log (2* (rhok+rho) )  — 

+log(x/2/i/alpha)*  ... 

(i-l)/2+log(abs( jA (i-1) *bessel (i-1, j*sqrt (2*i*alpha*x) 

/(rhok+rho)))) ; 

g=exp (g) *delta ; 

f i=f ind(isnan(g) ) ; 

g(fi)=zeros(l, length (fi) ) ; 

if  all(g==0)==l, 

gmin=l; 

else, 

[xx,gmax]=max(g) ; 

g(l:gmax)=ones(l,gmax) *xx; 

[xx,gmin]=min( (g-1.0e-6) .  A2) ; 

end 

wlimx=gmin*delta; 

end 

rhok=rhok2 ; 

delta=limx/100 ; 

x=0 . 5*delta : delta : 1 imx-0 . 5*delta ; 

g=- (x+2* (L-i) *alpha)/2/ (rhok+rho) -log (2* (rhok+rho) ) 

+log(x/2/(L-i)/alpha)*    ... 
(L-i-l)/2+log(abs(jA(L-i-l)*bessel(L-i-l, j*sqrt (2* (L-i) 

*alpha*x)/    . . . 
(rhok+rho) ) ) ) ; 
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g=exp (g) *delta ; 

f i=find(isnan(g) ) ; 

g (fi)=zeros(l, length (fi) )  ; 

if  all(g==0)==l, 

gmin=l; 

else, 

[xx,gmax]=max(g) ; 

g ( 1 : gmax ) =ones ( 1 ,  gmax ) *xx ; 

[xx,gmin]=min( (g-1.0e-6) .  A2) ; 

end 

zlimx=gmin*delta ; 

check  the  sequence  limit  again 
while  gmin  <=2, 
delta=zlimx/lOO; 

x=0 . 5*delta : delta : zlimx-0 . 5*delta ; 

g=-(x+2*(L-i)*alpha)/2/(rhok+rho)-log(2*(rhok+rho) ) 
+log(x/2/(L-i) /alpha)*  ... 

(L-i-l)/2+log(abs(jA(L-i-l)*bessel(L-i-l, j*sqrt (2* (L-i) 
*alpha*x) / (rhok+rho) ) ) ) ; 
g=exp(g) *delta; 
f i=find(isnan(g) ) ; 
g(fi)=zeros(l/ length (fi) ) ; 
if  all(g==0)==l, 
gmin=l; 
else, 

[xx/gmax]=max(g) ; 
g(l:gmax)=ones (1, gmax) *xx; 
[xx/gmin]=min( (g-1.0e-6) . A2)  ; 
end 

zlimx=gmin*delta ; 
end 

produce  the  sequency  of  jammed  hop  of  signal  branch 
del ta=wl imx/1 im ; 

x=0.49*delta:delta:wlimx-0. 51*delta; 

gl=-(x+2*i*alpha)/2/(rhokl+rho)-log(2*(rhokl+rho) )  ... 
+log (x/2/i/alpha) *  . . . 

(i-l)/2+log(abs(jA(i-l)*bessel(i-l, j*sqrt(2*i*alpha*x) 
/(rhokl+rho) ) ) ) ; 
gl=exp(gl) *delta; 
f i=f ind(isnan(gl) ) ; 
gl(fi)=zeros(l, length (fi) ) ; 
if  all(gl==0)==l/ 
gl=zeros(l, lim) ; 
gl(l)=l; 
else, 

gl=gl/sum(gl) ; 
end 

produce  the  sequence  of  non jammed  hop  of  signal  branch 
if  wlimx  >=  100*zlimx 
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gs=[gl  zpad] ; 

else 

g2=-(x+2*(L-i)*alpha)/2/(rhok2+rho)-log(2*(rhok2+rho) )  ... 

+log(x/2/(L-i)  

/alpha) * (L-i-l)/2+log(abs (jA (L-i-1) *bessel (L-i-1, 

j*sqrt(2*(L-i)  ... 

*alpha*x) / (rhok2+rho) ) ) )  ; 

g2=exp (g2 ) *delta ; 

fi=find(isnan(g2) ) ; 

g2 ( fi)=zeros(l, length (fi) ) ; 

if  all(g2==0)==l, 

g2=zeros (1,11m) ; 

g2(l)=l; 

else, 

g2=g2/sum(g2)  ; 

end 

%  produce  the  sequence  of  signal  branch 
gs  =  abs(ifft(fft(gl,n) .*fft(g2,n) ,n) ) ; 
gs  =  gs/sum(gs) ; 
end 

%  produce  the  sequence  of  jammed  hop  of  noise  branch 

pl=-log(2*rhokl)*i+log(x)*(i-l) -log (fact (i-1) ) -x/2/rhokl; 

pl=exp(pl) *delta; 

f i=f ind(isnan(pl) ) ; 

pi (fi)=zeros(l/ length (fi) ) ; 

if  all(pl==0)==l, 

pl=zeros(l/ lim) ; 

pl(l)=l; 

else, 

pl=pl/sum(pl) ; 

end 
%  produce  the  sequence  of  non jammed  hop  of  noise  branch 

p2=-log (2*rhok2) * (L-i) +log (x) * (L-i-1) -log (fact (L-i-1) ) -x/2/rhok2 ; 

p2=exp (p2 ) *delta ; 

fi=find(isnan(p2) ) ; 

p2 (fi)=zeros(l,length(fi) )  ; 

if  all(p2==0)==l/ 

p2=zeros(l/lim) ; 

p2(l)=l; 

else, 

p2=p2/sum(p2)  ; 

end 
%  produce  the  sequence  of  noise  branch 

ps  =  absfifftfffttp^n)  .*fft(p2,n)  ,n) )  ; 

ps  =  cumsum(ps) ; 

ps  =  ps/max(ps) ; 

end 
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%  compute  the  probability  of  bit  error  using  sequence 
tdata( j5+l)=tdata( j5+l) +bino(L, i) *gamaAi* (1-gama) A (L-i) 
*M/2*(l-sum(gs.*ps.^(M-l)) ) ; 
end,  end,  end 

data= [ data ; td ] ; 

save  datafile  data; 

end , end , end , end , end 
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APPENDIX  B 


%  This  is  the  program  to  produce  the  probability  of  bit  error  of 
%  Self-Normalization  FH/MFSK  Receiver 
%  Lee,  Kang  Yeun  (Major,  ROKAF) 

%  three  special  routines 

%  fact  :  factorial  routine 

%  bino  :  binomial  coefficient  routine 

%  ereval  :  error  probability  computation  routine  using  fft 

%  3  signal-to-noise  ratio (13 .35db,  16db,  20db) 

SNR=[21.62718524D0  39 . 81071706D0  100. 0D0] ; 
%  3  direct-to-diffuse  signal  power  ratio  (0.01,  1,  10) 

DDR=[0.01D0  1.0D0  10.0D0]; 
%  4  jamming  ratio 

JRATIOl=[1.0D0  0.2D0  0.1D0  0.05D0  0.01D0]; 

JRATIO2=[1.0D0  0.2D0  0.05D0  0.01D0  0.001D0]; 

%  Bandwidth-to-Symbol  Rate  Ratio  (BW/Rs  :  1   2  ) 
mu=1.0D0; 

%  M-ary 

M=16.0D0; 
%  Multichannel  factor  w=logM/log2 

omega=log(M)/log(2.0D0) ; 

%  Bit  SNR   13.35  db 
nu=SNR(l) ; 

count =1 ; 

%  Hopping  number  (1  to  6) 
for  Jl=l:6 
L=J1; 

%  Direct-to-Diffuse  signal  ratio  (0.01   1   10  ) 
for  J3=l:3 
psi=DDR(J3) ; 

%  Direct  SNR  of  NOJAMMING  CASE 

rho=omega*nu*psi/mu/(psi+1.0D0)/L; 
%  Diffuse  SNR  of  NOJAMMING  CASE 

xi=omega*nu/mu/ (psi+1. 0D0)/L; 

%  produce  the  sequence  of  no jamming  case 
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Z=-0.999D0:0.002D0:0.999D0; 

NOJAM=exp(rho*(z-1.0D0) ./ (2 . ODO+xi* (1 . ODO-z) ) ) 

*(1.0DO+z) .A (M-2. ODO) . . . 

./ (2. ODO+xi* (1. ODO-z) ) .*( (2 . ODO* (1. ODO+xi) ) ./ 

(xi*(1.0DO-z)+2.0DO 

.*(1.0DO+rho*(1.0DO+z) ./( (2 . ODO+xi* (1 . ODO-z) ) 

*(1. ODO+xi) ) )+M-2.0D0) ; 


z=0.001D0:0.002D0:0.999D0; 

B=rho*z./( (1.0DO+xi*(l. ODO-z) )*(1. ODO+xi) ) ; 

E=exp(rho*(z-1.0DO) ./ (1. ODO+xi* (1. ODO-z) ) ) ; 

TT=zeros(M-l,500) ; 

for  N=0:M-2 

TT(N+1, : )=fact(M-2.0D0)/fact(M-2.0D0-N)/fact(N) A2*B. AN; 

end 

Tl=sum(TT) .*((1.0DO-z)*(1.0DO+xi) ) .A(M-2) .*z.  ... 

*  (M-2. ODO) ./(l.ODO+z) . . . 

./(I. ODO+xi* (1. ODO-z) ) .A(M-1) .  *E; 

TT=zeros (M-2 , 500) ; 

for  N=0:M-3 

TTl=zeros(M-l-N,500)  ; 

for  IQ=0:M-2-N 

TT1 (IQ+1 , : ) =f act (M-2 . ODO-N) /f act (M-2 . ODO-N-IQ) /f act (IQ) A2*B. AIQ ; 

end 

TEMP=2.0D0*( (1. ODO-Z) *(1. ODO+xi) ) .A(M-3-N) . *(1. ODO+xi) .*z.A2 

.* (l.ODO+z) .A(N-1) . * (M-2. ODO-N) ; 
TEMP=TEMP./(1. ODO+xi* (1. ODO-z) ) . A (M-l . ODO-N) . *E. *sum(TTl) ; 
TTl=zeros(M-2-N,500) ; 
for  IQ=0:M-3-N 

TT1 (IQ+1 , : ) =f act (M-3 . ODO-N) /fact (M-3 . ODO-N-IQ) /fact (IQ) A2*B. AIQ; 
end 
TEMP1=( ( 1. ODO-z) *(1. ODO+xi) ) .A(M-3-N) . *z . * (1 . ODO+z)  . . . 

.AN.* (M-2. ODO) ./(l. ODO+xi* (1. ODO-Z) ) .A(M-2-N) .*E; 
if  IQ  ==  0, 

TT (N+l , : ) =TEMP+TEMP1 . *TT1 ; 
else 

TT(N+1, :)=TEMP+TEMPl.*sum(TTl) ; 
end 
end 
Tl=Tl+sum(TT) ; 

%  sequence  of  nojamming  case 

NOJAM (501: 1000) =NOJAM( 501: 1000) -Tl; 

Ps=zeros (5,41) ; 

%  jamming  ratio 
for  J4=l:5 
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if  J3==3 

gamma=JRATI02 (J4) ; 
else 

gamma=JRATI01(J4) ; 
end 

%  Signal-to-Jamming  noise  ratio  (  0  db  -  40  db  ) 

for  J5=0:40 
iota=10.0D0A(0.1D0*J5) ; 

%  Direct  SNR  of  JAMMING  case 

rho=omega*nu*iota*psi*gamma/mu/ (psi+1. 0D0)/ (iota*gamma+nu)/L; 
%  Diffuse  SNR  of  JAMMING  case 

xi=omega*nu*iota*gamma/mu/ (psi+1. 0D0)/ (iota*gamma+nu)/L; 

%  produce  the  sequence  when  the  signal  is  jammed  by  noise  jamming 

z=-0.999D0:0.002D0:0.999D0; 

JAM=exp(rho*(z-1.0D0) ./ (2 . ODO+xi* (l.ODO-z) ) ) .  ... 
*(1.0D0+z) .A(M-2.0D0) . . . 

./(2.0D0+xi*(1.0D0-z)) .  ... 
*(  (2. 0D0*(1. ODO+xi) ) ./ (xi*( l.ODO-z) +2. 0D0) . . . 
,*(1.0D0+rho*(1.0D0+z) ./( (2. ODO+xi* (1. ODO-z)  )  ... 
*(1. ODO+xi) ) )+M-2.0D0) ; 

Z=0.001D0;0.002D0:0.999D0; 

B=rho*z./( (1.0D0+xi*(l. ODO-z) )*(1. ODO+xi) ) ; 

E=exp(rho*(z-1.0D0) ./ (1. ODO+xi* (1. ODO-z) ) ) ; 

TT=zeros(M-l,500) ; 

for  N=0:M-2 

TT(N+1, : )=f act (M-2.0D0) /fact (M-2.0D0-N) /fact (N) A2*B.AN; 

end 

Tl=sum(TT) .*( (l.ODO-z) *(1. ODO+xi) ) .A(M-2) .  ... 

*z.*(M-2.0D0) ./(l.ODO+z) . .  . 

./(I. ODO+xi* (l.ODO-z)) .A(M-1) .*E; 

TT=zeros (M-2 , 500) ; 

for  N=0:M-3 

TTl=zeros(M-l-N/500) ; 

for  IQ=0:M-2-N 

TT1 (IQ+1 , : ) =f act (M-2 . 0D0-N)/f act (M-2 . ODO-N-IQ) /f act (IQ) A2*B. AIQ; 

end 

TEMP=2.0D0*( (l.ODO-z) *(1. ODO+xi) ) .A(M-3-N) . *(1. ODO+xi) .*z.A2  ... 

.* (l.ODO+z) .A(N-1) .*(M-2.0D0-N) ; 
TEMP=TEMP./(1. ODO+xi* (l.ODO-z) ) . A (M-l. ODO-N) . *E. *sum(TTl) ; 
TTl=zeros(M-2-N,500) ; 
for  IQ=0:M-3-N 
TTltlQ+l, : )=f act (M-3. ODO-N) /f act (M-3. ODO-N-IQ) /f act (IQ)A2*B.AIQ; 
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end 

TEMPl=((1.0D0-z)*(1.0D0+xi) ) . A(M-3-N) . *z . * (1 . ODO+z)  . . . 

.AN.*(M-2.0D0) ./(l-0D0+xi*(1.0D0-z) ) . *(M-2-N) .  *E; 
if  IQ  ==  0, 

TT (N+l , : ) =TEMP+TEMP1 . *TT1 ; 
else 

TT(N+1, :)=TEMP+TEMPl.*sum(TTl) ; 
end 
end 
Tl=Tl+sum(TT) ; 

%  the  sequence  of  jamming  case 
JAM (501: 1000) =JAM( 501: 1000) -Tl ; 

%  'bino'  and  'ereval'  are  the  functions  for  binomial  and  error 
%   evaluation. 

%   Ps(i)  <=  (M-l)  *  Pr(  z  <  0  ) 

%   Pb  =  (M/2)/(M-l)*Ps 

%   Pb  <=  (M/2)  *  Pr(  z  <  0  ) 

%  probability  of  symbol  error 

for  i=0:L 

Ps(J4,J5+l)=Ps(J4,J5+l)+bino(L, i) *gammaAi* (1-gamma) A (L-i)  . . . 

*ereval(NOJAM,L-i,JAM,i) ; 
end 

%  probability  of  bit  error 
Ps(J4/J5+l)=(M/2.0D0)*Ps(J4/J5+l) ; 

end 
end 

%  Save  the  Results  and  Plot 

save  data  Ps 

end 
end 
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